Thermal neutron captures on d and 3 He 
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We report on a study of the nd and n 3 He radiative captures at thermal neutron energies, using 
wave functions obtained from either chiral or conventional two- and three-nucleon realistic potentials 
with the hyperspherical harmonics method, and electromagnetic currents derived in chiral effective 
field theory up to one loop. The predicted nd and n 3 He cross sections are in good agreement with 
data, but exhibit a significant dependence on the input Hamiltonian. A comparison is also made 
between these and new results for the nd and n 3 He cross sections obtained in the conventional 
framework for both potentials and currents. 
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The nd and n 3 He radiative capture reactions at ther- 
mal neutron energy are very interesting, in that the mag- 
netic dipole (Ml) transitions connecting the continuum 
states to the hydrogen and helium bound states are inhib- 
ited at the one-body level. Hence, most of the calculated 
cross sections (80-90% in the case of n 3 He) results from 
contributions of many-body components in the electro- 
magnetic current operator [l| . Thus these processes pro- 
vide a crucial testing ground for models describing these 
many-body operators and, indirectly, the nuclear poten- 
tials from which the ground- and scattering-state wave 
functions are derived. 

Over the past two decades, chiral effective field theory 
(xEFT), originally proposed by Weinberg in a series of 
papers in the early nineties Q , has blossomed into a very 
active field of research. The chiral symmetry exhibited 
by quantum chromodynamics (QCD) severely restricts 
the form of the interactions of pions among themselves 
and with other particles. In particular, the pion couples 
to baryons, such as nucleons and A-isobars, by powers of 
its momentum Q, and the Lagrangian describing these 
interactions can be expanded in powers of Q/A x , where 
A x ~ 1 GeV specifies the chiral-symmetry breaking scale. 
As a result, classes of Lagrangians emerge, each charac- 
terized by a given power of Q/A x and each involving a 
certain number of unknown coefficients, so called low- 
energy constants (LEC's), which are then determined by 
fits to experimental data (see, for example, the review 
papers 0| and and references therein). Thus, xEFT 
provides, on the one hand, a direct connection between 
QCD and its symmetries, in particular chiral symmetry, 
and the strong and electroweak interactions in nuclei, 
and, on the other hand, a practical calculational scheme 
susceptible, in principle, of systematic improvement. In 
this sense, it can be justifiably argued to have put low- 
energy few-nucleon physics on a more fundamental basis. 

Concurrent with these conceptual developments have 
been the acquisition and refinement of accurate meth- 



ods for solving the A=3 and 4 Schrodingcr equation (see 
Ref. Q for a review). In this respect, it is worthwhile 
noting that the A=A scattering problem has proven to 
be especially challenging for two reasons. The first is 
its coupled-channel nature: even at vanishing energies 
for the incident neutron, the elastic n- 3 He and charge- 
exchange p- 3 H channels are both open, and need to be 
accounted for. The second reason lies in the peculiari- 
ties of the 4 He spectrum, in particular the presence of 
resonant states between the p- 3 H and n- 3 Hc thresholds, 
which make it hard to obtain numerically converged so- 
lutions. Indeed, it is only very recently that both these 
capabilities have been fully realized @, 0| • In the present 
work, the 3- and 4-body problems are solved with the 
hyperspherical-harmonics (HH) technique (see Ref. Q for 
a review). 

The developments outlined above make it possible to 
examine the question of whether available experimen- 
tal data on these delicate processes — the nd and n 3 He 
captures — are well reproduced by theory. The present 
letter reports on such an effort by presenting results ob- 
tained both in xEFT as well as in the conventional frame- 
work based (essentially) on a meson-exchange model of 
potentials and electromagnetic current operators. This 
approach, while more phenomcnological than xEFT, has 
a broader range of applicability, and accounts satisfacto- 
rily for a wide variety of nuclear properties and reactions 
up to energies, in some cases, beyond the pion production 
threshold, see Ref. [H[ for a review. In particular, it re- 
produces well observed magnetic properties of A = 2 and 
3 nuclei, including moments and form factors, as well as 
the np radiative capture, see Marcucci et al. (2005) [l|. 

The model for the nuclear electromagnetic current in 
%EFT up to one loop was derived originally by Park 
et al. using covariant perturbation theory. In the 
last couple of years, two independent derivations, based 
on time-ordered perturbation theory (TOPT) , have ap- 
peared in the literature, one by Kolling et al. [101 ] and the 
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other by some of the present authors [ll| . There are tech- 
nical differences in the implementation of TOPT, which 
relate to the treatment of reducible diagrams and are doc- 
umented in considerable detail in the abo ve p apers. How- 
ever, the resulting expressions in Refs. and for 
the two-pion-exchange currents (the only ones considered 
by the authors of Ref. [lo|) are in agreement with each 
other, but differ from those of Ref. [9j, in particular in 
the isospin structure of the Ml operator associated with 
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the one-loop corrections — see Pastore et al. (2009) [ll| 
for a comparison and analysis of these differences. 

Explicit expressions for the %EFT curren t s U p t G 
one loop, and associated Ml operators, arc listed in 
Refs. [11] . Here we summarize succinctly their main fea- 
tures. The leading-order (LO) term results from the 
coupling of the external photon field to the individ- 
ual nucleons, and is counted as eQ~ 2 (e is the electric 
charge). The NLO term (of order eQ^ 1 ) involves seag- 
ull and in-flight contributions associated with one-pion 
exchange, and the N 2 LO term (of order e Q°) represents 
the (Q/tojv) 2 rclativistic correction to the LO one-body 
current (m^r denotes the nucleon mass). 



(a) (b) (c) (d) (e) (f) (g) (h) (i) (j) (k) 

FIG. 1. Diagrams illustrating two-body currents at N 3 LO. Nucleons, pions, and photons are denoted by solid, dashed, and 
wavy lines, respectively. Only one among the possible time orderings is shown for diagrams (a)-(j). 



At N 3 LO (e Q) we distinguish three classes of 



terms ll| : i) two-pion exchange currents at one loop, 
illustrated by diagrams (a)-(i) in Fig. [1] ii) a tree-level 
one-pion exchange current involving the standard irNN 
vertex on one nucleon, and a jirNN vertex of order e Q 2 
on the other nucleon, illustrated by diagram (j), and hi) 
currents generated by minimal substitution in the four- 
nucleon contact interactions involving two gradients of 
the nucleons' fields as well as by non-minimal couplings, 
collectively represented by diagram (k). A fourth class 
consisting of (Q/uin) 2 relativistic corrections (RC's) to 
the NLO currents is neglected. However, RC's are not 
consistently treated in available chiral potentials, such as 
those employed below. For example, the RC's in the two- 
nucleon potential, implied by Poincare covariance and 
just derived in Ref. [l2| in an EFT context, have been 
omitted so far in A=3 and 4 calculations, even though 
their contribution is expected to be comparable to that 
of the three- nucleon potential 12 1. 



The loop corrections in panels (a)-(i) involve the 
pion mass, the nucleon axial coupling constant (7^4 =1.29 
(from the Goldbcrgcr-Trciman relation relating it to the 
wNN coupling constant), and the pion decay amplitude 
F OT =184.8 MeV. The LEC's entering panel (i) and the 
minimal currents in panel (k) have been determined by 
fits to the np S- and P-wave phase shifts up to 100 MeV 
laboratory energies [ll|. We refer below to these con- 
strained terms collectively as N 3 LO(S-L). There are five 
additional unknown LEC's: dg, dg, and c^i m panel (j), 
and C' 15 and C( 6 in the non-minimal currents of panel 
(k). We denote these terms as N 3 LO(LECs) in the fol- 
lowing. In a resonance saturation picture, the c?g and d^i 
(dg) LEC's can be related to the combination of coupling 
constants and masses entering the isovector (isoscalar) 



N-A excitation and uittj (pirj) transition currents 11 
Indeed, this connection is exploited in a series of calcula- 
tions, based on the Ml operators derived in Ref. 0], of 
the np, nd, and n 3 Hc radiative captures, and magnetic 
moments of A~2 and 3 nuclei 0, [l3|. Here, however, we 
adopt a different strategy, as discussed below. Lastly, we 
observe that at N 3 LO there are no three-body currents 
in the formalism of Ref. 11], which retains irreducible 
and recoil-corrected reducible diagrams. 

The %EFT Ml operators have power- law behavior 
for large relative momenta fc's, and need to be reg- 
ularized, before they can be inserted between nuclear 
wave functions. Following common practice, we imple- 
ment this regularization by means of a cutoff C\(k) = 
cxp(-/c 4 /A 4 ), with A in the range (500-700) MeV, and 
constrain the LEC's entering the N 3 LO Ml operators of 
panels (j) and (k) in Fig. [T] to reproduce a set of observ- 
ables for any given A in this range. This same renor- 
malization procedure is adopted in the currently most 
advanced analyses of nuclear potentials, for example in 
Ref. [3 (see also Sec. II. C of Ref. and references 
therein, for further discussion of this issue). 

These operators are used in the present work to study 
the magnetic moments of the deuteron and trinucleons, 
and the np, nd, and n 3 He radiative captures at ther- 
mal neutron energies. The calculations are carried out 
by evaluating their matrix elements between wave func- 
tions obtained from either conventional or chiral (re- 
alistic) potentials with the variational HH method 
We consider the Argonne V\g [l5| (AV18) and chiral 
N 3 LO [3] (N3LO) two-nucleon potentials in combina- 
tion with the Urbana-IX (UIX) and chiral N 2 LO [13] 
(N2LO) three-nucleon potentials. The AV18/UIX and 
N3LO/N2LO Hamiltonians provide a very good descrip- 



tion of three- and four-nucleon bound and scattering 
state properties, including binding energies, radii, and 
scattering lengths 0, . 

We now turn our attention to the determination of the 
LEC's dg, dg, d' 2l , C 15 , and C 16 . In principle, the d- could 
be fitted to pion photoproduction data on a single nu- 
cleon or, as mentioned already, related to hadronic cou- 
pling constants by resonance saturation arguments. Both 
procedures have drawbacks. While the former achieves 
consistency with the single-nucleon sector, it neverthe- 
less relies on single-nucleon data involving photon en- 
ergies much higher than those relevant to the thresh- 
old processes under consideration and real (in contrast 
to virtual) pions (some of these same issues in the con- 
text of three-nucleon potentials have been investigated in 
Ref. [IH ). The second procedure is questionable because 
of poor knowledge of some of the hadronic couplings, such 
as QujN n and g P NN- 

Here, we assume d^/d^ = 1/4 as suggested by A- 
dominance, and rely on nuclear data to constrain the re- 
maining four LEC's. The values obtained by reproducing 
the experimental np cross section and magnetic moments 
of the deuteron and trinucleons are listed in Table HI 
Note that the adimensional values reported there are in 
units of powers of A, i.e., we have defined rfg = e?f/A 2 , 
C' 15 = df/A 4 , d£ x = </A 2 , and C' w = 47 A 4 and the 
superscripts S and V denote the isoscalar and isovector 
content of the associated operators. 

TABLE I. Adimensional values of the LEC's corresponding 
to cutoff parameters A in the range 500-700 MeV obtained 
for the AV18/UIX (N3LO/N2LO) Hamiltonian. See text for 
explanation. 
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-8.85 (-0.225) 
-2.90 (9.20) 
6.64 (20.4) 


-3.18 (-2.38) 
-7.10 (-5.30) 
-13.2 (-9.83) 


5.18 (5.82) 
6.55 (6.85) 
8.24 (8.27) 


-11.3 (-11.4) 
-12.9 (-23.3) 
-1.70 (-46.2) 



In Fig. [5] we show results obtained by including cu- 
mulatively the contributions at LO, NLO, N 2 LO, and 
N 3 LO(S-L) for the deuteron (fid) and 3 Hc/ 3 H isoscalar 
(/i S ) magnetic moments (left panels), and for the np ra- 
diative capture cross section (a^ ip ) at thermal energies 
and 3 He/ 3 H isovector (/jY) magnetic moment (right pan- 
els). The NLO and N 3 LO(S-L) Ml operators are purely 
isovector, and hence do not contribute to fid and fi s . The 
band represents the spread in the calculated values corre- 
sponding to the two Hamiltonian models considered here 
(AV18/UIX and N3LO/N2LO). The sensitivity to short- 
range mechanisms, encoded in the cutoff C\(k) and in the 
rather different short-range behaviors of the adopted po- 
tentials, remains quite weak for all these observables. Of 
course, taking into account the N 3 LO(LECs) contribu- 
tion with the LEC values listed in Table U reproduces the 
experimental data represented by the black band (to ac- 
commodate errors, but these are negligible in the present 
case). The contributions at LO and NLO have the same 
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FIG. 2. Results for the deuteron and trinucleon isoscalar and 
isovector magnetic moments, and np radiative capture, ob- 
tained by including cumulatively the LO, NLO, N 2 LO, and 
N 3 LO(S-L) contributions. See text for discussion. 

sign, while those at N 2 LO and N 3 LO(S-L) have each op- 
posite sign, and tend to increase the difference between 
theory and experiment. 

Having fully constrained the x^FT Ml operator up 
to N 3 LO, we are now in a position to present predic- 
tions, shown in Fig. [31 for the nd and n 3 He radiative 
capture cross sections, denoted as a^ d and cr^ 3He , and 
the photon circular polarization parameter R c resulting 
from the capture of polarized neutrons on deuterons. The 
experimental data (black bands) are from Ref. [l9j for 
nd and Ref. [2(| for n 3 He. Results obtained with the 
complete N 3 LO xEFT operator are shown by the or- 
ange band labeled N 3 LO(LECs): those corresponding 
to the AV18/UIX (N3LO/N2LO) model delimit the up- 
per (lower) end of the band in the case of nd, and its 
lower (upper) end in the case of n 3 He. Their sensitiv- 
ity to the cutoff, within a given model, is negligible for 
nd and at the 5-10% level for n 3 Hc. The AV18/UIX 
and N3LO/N2LO results are within ~ 2% of the nd 
experimental cross section. However, at A=600 MeV, 
for example, the experimental cr^3 He is well reproduced 
in the N3LO/N2LO calculation, but underpredicted by 
~ 15% in the AV18/UIX. As expected, these processes 
are strongly suppressed at LO: the calculated <j^ d (LO) 
and c^Hc^O) are l ess than half and a factor of five 
smaller than the measured values. In the case of n 3 Hc, 
the matrix element at NLO is of opposite sign and twice 
as large (in magnitude) compared to that at LO, hence 
cr^3 He (LO) and <7^ 3He (LO + NLO) are about the same, 
as seen in Fig. [3j For nd, however, the LO and NLO 
contributions interfere constructively. For both nd and 
n 3 He, the N 2 LO and N 3 LO(S-L) corrections exhibit the 
same pattern discussed in connection with Fig. [2l The 
N 3 LO(LECs) contributions are large, and essential for 
bringing theory into good agreement with experiment. 

In Fig. [3J we also show results obtained in the con- 
ventional framework, referred to as the standard nuclear 
physics approach (SNPA), with the AV18/UIX Hamil- 
tonian model. The electromagnetic current operator in- 
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FIG. 3. Results for a~^ d (left top panel), o"^3 He (right 
top panel), and R c (left bottom panel), obtained by in- 
cluding cumulatively the LO, NLO, N 2 LO, N 3 LO(S-L), and 
N 3 LO(LECs) contributions. Also shown are predictions ob- 
tained in the standard nuclear physics approach (squares la- 
beled SNPA and SNPA*). See text for discussion. 



eludes the one-body term — the same as the x^FT LO op- 
erator discussed earlier — as well as two- and three-body 
terms, constructed from the two- and three-nucleon po- 
tentials (AV18 and UIX, respectively) so as to satisfy 
exactly current conservation with them, see Marcucci et 
al. (2005) [J. In the figure, the squares labeled SNPA* 
represent the results obtained by retaining in addition 
the relativistic corrections to the one-body current (i.e., 
the xEFT N 2 LO operator). These corrections had been 
neglected in all previous studies of these processes [l[ . In 
fact, their contributions are found to be numerically sig- 
nificant and, at least for the case of nd capture, bring the 
present SNPA* results within ~ 4% and 2%, respectively, 
of the experimental data and %EFT predictions (based 
on the AV18/UIX model). However, it should be em- 
phasized that the SNPA (and SNPA*) currents contain 
no free parameters — i.e., they are not constrained to fit 
any photonuclcar data, in contrast to the %EFT currents. 
From this perspective, the achieved level of agreement be- 
tween SNPA* and data should be viewed as satisfactory, 
especially when considering, in the ^EFT context, the 
large role played by the N 3 LO(LECs) currents. 

We conclude by remarking that the convergence of the 
chiral expansion is problematic for these processes. The 
LO is unnaturally small, since the associated operator 
cannot connect the dominant S-states in the hydrogen 
and helium bound states (in contrast to np capture, for 
example) This leads to an enhancement of the NLO, 
which, however, in the case of n 3 He is offset by the de- 



structive (and accidental) interference between it and the 
LO contribution. It appears that at N 4 LO no additional 
LEC's enter, Park et al. (2000) 0. Thus, inclusion of 
the N 4 LO currents would have to be followed by a "reseat- 
ing" of the LEC's in Table [I] in order to reproduce (as 
before at N 3 LO) the experimental values of fid, f- S ' V , and 
a^[ p . The resulting predictions for a^ d and T^ 3He would 
presumably be close to those obtained here at N 3 LO. 

It is likely that explicit inclusion of A degrees of free- 
dom would significantly improve the convergence pattern, 
particularly in view of the relevance of the pion-exchange 
current of panel (j) in Fig.Q] — in such a theory, this oper- 
ator would be promoted to N 2 LO [11]. We plan to pursue 
vigorously this line of research in the future. 

The work of R.S. is supported by the U.S. Department 
of Energy, Office of Nuclear Physics under contract DE- 
AC05-06OR23177. 
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